KNK 962 



Two-Body Dirac Equation and Its Wave Function 

at the Origin * 

Hitoshi Ito 

Department of Physics, Faculty of Science and Engineering 
Kinki University, Higashi- Osaka, 577 Japan 

December 18, 1996 



Abstract 

A relativistic equation is deduced for the bound state of two particles, by assuming a 
proper boundary condition for the propagation of the negative-energy states. It reduces to 
the (one-body)Dirac equation in the infinite hmit of one of the constituent mass. It also 
has the symmetries to assure the existence of the anti-bound-state with the same mass. 
The interaction kernel(pseudo-potential) is systematically constructed by diagonalizing the 
Hamiltonian of the background field theory, by which the retardation effects are included 
in the interaction. Its wave function at the origin(WFO) behaves properly in a manner 
suggested by the covariant field theory. The equation can well be applied in the heavy 
quark effective theory. 
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1 Introduction 



One of the unsatisfactory nature of the Bethe-Salpeter equation for the bound state 
is that it does not reduce to the Dirac equation in the infinite hmit of the one of the 
constituent mass, even when the interaction is assumed to be instantaneous|]l[. We 
have to sum up all the crossed diagrams to recover the Dirac equation^], which is 
impossible for the finite masses. The reducibility to the Dirac equation is essential 
in applications to the heavy flavored quarkonia or the other atom-like systems. 

Historically, the relativistic single-time equation for the two-body system pre- 
ceded the BS equation. Soon after the discovery of the Dirac equation, Breit pro- 
posed the equation of the form|Q 

{H,{p,) + H2{p,) + V}i^ = Ei^, (1) 
where Hi is the Dirac Hamiltonian 

Hi{pi) = OLi -Pi + rriif^i 

and \^ is a local potential. The Breit equation reduces to the Dirac equation in 
the limit mentioned above but does not have the "i?-parity symmetry", by which 
we mean there is symmetric negative eigenvalue E for each positive one, which is 
interpreted as the bound state of the antiparticles. 

E'-parity symmetry is a result of the PCT invariance. The Dirac Hamiltonian 
is odd under the PCT transformation and the interaction Hamiltonian considered 
below is invariant under it. We therefore see that the i?-parity symmetry in the 
instantaneous BS equation is assured by a projection factor 

A++ - A- (2) 
in front of the potential V, which consists of the energy-projection operators 

Ae„(Pi,P2) = A^(Pi)Aj(P2), e,r] = + oT-, (3) 

where 

Alip,) = {EM+eHiip,)}/2Eiip,), 

E^{p^) = {pl + mlfl\ 

We must introduce a similar factor in any attempt at the construction of an improved 
two-body equation. 

The factor comes from the Stiickelberg-Feynman boundary condition for the 
propagation of the negative-energy state [Q. We will construct the equation for the 
unequal-mass constituents by imposing similar boundary condition. Before pre- 
senting it, we restrict the framework of consideration. For definiteness, we assume 
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Abelian gauge field interacting with the Dirac particles. We also work in the rest 
frame of the bound system(P = {E,0)), since we are looking for the non-covariant 
approximation of the low-energy dynamics, p and x, in the following, are the relative 
momentum and coordinate, respectively, in this frame. 



2 The equation for unequal-mass constituents 

In this section, we deduce a new single-time equation(Two-body Dirac equation) for 
the unequal- mass constituents, for which we assume that the mass M(mi) is larger 
than m(m2). 



2.1 Two-body Dirac equation 

Let us first consider the 2-body propagator for which we impose the boundary condi- 
tion that the negative-energy states propagate backward in time. If we keep individ- 
uality of the constituents and use the usual Feynman propagator in the instantaneous 
BS equation, the equation (p modified by the factor results. However, we can 
choose the other possibility in which we incorporate the idea that the bound two 
bodies should be treated as a quantum-mechanical unity.Q Since M is larger than m, 
the free part of the Hamiltonian has the same sign as that of the particle 1 in the CM 
system. Let us then modify the boundary condition as follows: A bound two-particle 
state propagates backward in time if their net energy is negative. By assuming it, we 
have the propagator 



where Ai -|- A2 = 1 and the limit 5 — > +0 is assumed. 

By assuming (^) we obtain the single-time equation. The pseudo-4-dimensional 
form of it, in the momentum space, is 

ij{p) = iSf{P,p) J V{p, q)^{q)dy {27r)\ (5) 

where V{p, q) represents the interaction, which does not depend on the relative 
energies but is not necessarily an instanteneous local potential. After integrating out 
the redundant degree of the freedom, (||) becomes 

^We can establish the concept of individuaUty in the quantum mechanics only through observa- 
tion, which brings about a subtle point to the bound system: To detect an individual one in the 
bound constituents, we need to separate them by applying the 3rd interaction, which inevitably 
destroys the original state. So there is no a priori reason to apply the free propagator individually 
to each constituent. Note that "the free propagator" itself for the bound state is merely a convention 
of approximation. 
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{H,{-p) + H2{p)+Y.^KnV}^ = E^, (6) 

er) 

which reduces to the Dirac equation in the infinite hmit of M. It is easy to see the 
£^-parity symmetry of this equation. 

If we would apply our equation to the scattering state, we shall have, from the 
time-dependent equation, the conserved probability density 

p{x, t) = tp{x, ty ^^ev4^{x, t), (7) 

er] 

which is in accord with the boundary condition that the negative- energy state prop- 
agates backward in time carrying the negative probability density]^. However, it 
does not necessarily provide the normalization condition for the bound state: For 
the scattering processes, the projected wave function A^^tp corresponds to the physi- 
cal state of the (free) particles with the positive or negative energy E. And the above 
interpretation of the probability current actually says that the state with the negative 
E is carrying the negative probability density. However, for the bound state with the 

positive eigenvalue E, A or A^^ip is merely a negative- energy component in the 

representation in which the energy of the free particle is diagonal. 

Taking the above consideration into account, we restore the probability interpre- 
tation of the wave function and normalize it by assuming the probability density 

p{x) = ^{x,ty^{x,t). (8) 
Observables except the Hamiltonian are self-adjoint under this metric: 

(0,6^) = (00,^). (9) 

The Hamiltonian is the operator ruling the time development of the system and its 
interaction part is modified by the factor J2er]^Aeri. Though it is not a self-adjoint 
operator, its eigenvalue is proved to be real if the inner product (pTS]) below exists. 

2.2 Green's function and the vertex equation 

We define the Green's function G for (|^) by the operator equation 

{E-H,-H2-J2 eAer]V}G = J2 ^A.„ (10) 

er] er] 

and 

G{E-H^-H2-VJ2 ^^er]} = E ^^er]. (H) 

er] er] 
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In the momentum representation, it can be written, by using the eigen-f unction Xn{p) 
of (Q), as 

G{P, P') = TTTB ^Xn{p)Xn{p'Y + continuum, (12) 

where En is an eigenvalue and is a normahzation factor defined by 

Nn = {Xn^^eKriXn)- (13) 
erj 

When one of the constituents(labeled with 2) belongs to the same class of the antipar- 
ticle of the other, there can be an annihilation process for which the unamputated- 
decay-vertex $ is given by 

$ = C-iG, (14) 

where 7 is the lowest vertex and C is the charge- conjugation matrix of the particle 
2. 

$ satisfies the vertex equation 
and the amputated vertex is 

r = 7 + C^V. (16) 

We can determine the renormalization constant Zi for the wave function at the 
origin (WFO) from (|T5|) and (0), if we need it|^. 

2.3 Interaction Hamiltonian 

We have, so far, not specified the interaction Hamiltonian (quasipotential). In this 
section, we investigate it for the one-(Abelian)gauge-boson exchange in the Coulomb 
gauge as an example. For the instantaneous part of the interaction, V is obvious. 
For the remaining part, we can specify the quasipotential in a clear way from the 
background field theory. We have already shown, for the Salpeter equation, that the 
quasipotential from the one-boson exchange is given through the diagonalization of 
Fukuda, Sawada, Taketani[0 and Okubopl (FSTO)0. We also apply this method 
to the present equation. 

We introduce the generalized Fock subspace of the free constituents, the bases 
of which are denoted by |e, ?7,p), where e and 77 are the signs of the energies of 
the particles 1 and 2 respectively.^ We then diagonalize the Hamiltonian in the 

^There was some error concerning the Fock space in Ref.|^. Namely, we employed the 
usual (positive-energy) Fock space and reinterpreted the matrix elements of the interaction Hamilto- 
nian including the negative-energy indices as the ones in this space according to the hole theory. It, 
however, brings the procedure into confusion, since we have revived the negative energy in Eq.(^). 
However, the error is only conceptual for the Salpeter equation which has only projection factors 
A_(__|_ and A The result of Ref.[^ is correct. 
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Schrodinger picture by using the FSTO method. The second-order boson-exchange 
potential in this subspace is given by 
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where q = p — P - The retardation effects are included in this equation. 



3 Some properties 

3.1 WFO of the ^Sq state 

Let us next study some basic feature of the equation. When it is applied to the system 
in which the pair annihilation of the constituents can occur, an important physical 
quantity is the wave function at the origin(WFO). For example, the decay amplitude 
of the pseudo-scalar Qq meson via a weak boson is proportional to the average WFO 
Tr{757o'?/'(0)}, where ip{0) is the charge conjugated (with respect to the particle 2(g)) 
WFO. We investigate, in Appendix, the asymptotic behavior of the momentum- 
space wave function by using the method given in Ref. 0. We assume instantaneous 
exchange of a gauge bosonQ. The average WFO thus obtained is finite. This result is 
consistent with consideration on the covariant field theory. We note that the average 
WFO becomes divergent in the limit of the one-body Dirac equation, for which we 
have the renormalization procedure [ITU . 



There are many "two-body Dirac equation" proposed. An interesting one from 
the point of view of the present paper is the one by Mandelzweig and Wallace |TT 



They intended to include the effects of the higher-order interaction (the crossed 
Feynman diagram) and got an equation which has the proper one-body limit and 
the ii^-parity symmetry. An important difference from our equation is in the average 
WFO considered above. It is divergent in their equation if the transverse part of the 
gauge-boson exchange is included | 



3.2 The crossed diagram 

By modifiing the boundary condition in (|^), we have included some parts of the 
crossed diagram. To estimate this effect, we take up the 4th-order scattering ampli- 
tude in the Coulomb potential for the mass-shell particles. We are interested in the 
contribution of the intermediate state with the A_(__ projection. Now, let us compare 

■'The analysis in the Appendix cannot be applied to the retarded interaction. 
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the two amplitudes restricted to this intermediate state: the 'ladder diagram' (/) with 
the propagator (|D and the crossed diagram(c) with the usual Feynman propagator. 
We see that they are different by the following factors: 



m 



for the 'ladder diagram' and 



(18) 



(P2 - 9)2 + m2 - a(2) . (p'2 - g) - 



(P2 ~ + rn'^i—A + a/p22 + m"^ + J {p'2 — qY + 



(19) 



for the crossed diagram, where ^"^^ P'2 initial and final momentum of 

the lighter particle respectively and q is the loop momentum to be integrated. The 
common element A is given by 



A = ^p,^ + M2 - ^{p, - qf + M\ 

We see that // and Ic have the same signature and the order of magnitude for mo- 
menta much smaller than M. 

One would intend to proceed to the next order approximation in the interaction 
V . He will calculate the 4th-order crossed diagram by assuming ordinary Feynman 
propagator. He will, further, have to add counter correction terms to compensate 
the part of the crossed diagram already included effectively in the iteration of the 
lowest kernel V , an example of which is given by the term with the factor J;. 



3.3 On application to the heavy quark effective theory 

One of the research fields in which we can utilize the two-body Dirac equation is the 
physics of the heavy flavored quarkonia. The recent trends in this field are led by 
the heavy quark effective theory (HQET). The Qq system is well described by using 
the two-body Dirac equation with some phenomenological (pseudo)potential. If we 
take the heavy quark limit (M 00) of the quark the equation itself becomes the 
one-body Dirac equation, which is the basic equation of HQET. There are, however, 
divergent(as M — >■ 00) portions in the WFO, some average of which determines the 
annihilation decay rate of a pseudoscalar quarkonium. The matrix element of the 
current for this decay is given by 

< 0|g7.75Q|Qg >~ C(^)tr{7,75^(0)}, (20) 

where ^(0) is the WFO of the bound state and we have separated the correction C{^) 
coming from the range of momenta /i ~ M in the integral over the loop momenta. 
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is usually determined from the perturbative loop corrections summed up by 
using the renormalization group equation |13l|14 . 



C(/^) = (^^r, rf= -6/(33 -2iV^). (21) 

We can calculate 1/M corrections to HQET by using the two-body Dirac equa- 
tion. 



3.4 Another approach to the heavy quark phenomenology 

An alternative way of investigating the heavy quark systems is due to direct appli- 
cation of the two-body Dirac equation. The binding interaction, in this formulation, 
consists of the Coulomb and the confining potentials. We further add the gluonic 
correction term to the potential, for which we should attach the high momentum 
cutoff [0. The correction factor C(/i) depends on this cutoff. If we take the cut- 
off of the QCD scale, we should include all the corrections except for the Coulomb 



contribution to the vertex correction. By subtracting it in the equation (pTD , we 



get d = 2/(33 — 2Nf) for the exponent [0] In the phenomenology of the cc and 
bb systems|P], we took fairly large value for this cutoff. If we choose the same one 
we only need the correction from the self-energy of the quarks. Then, the exponent 
becomes d = -2/(33 - 2A^/). 



3.5 On the equal-mass limit 

The unequal-mass equation (^) is well applied to the system with M ^ m. For 
M c:^ m, the reason justifiing (^ becomes obscure and we will have two different 
equations in the limit M = m. Though we concern ourselves in the case M > m, it 
is meaningful to investigate the degree of ambiguity near the equal mass limit. For 
M = m, the projection factor in front of the interaction term of (^ includes a part 
which violates the exchange symmetry: It is shown that 

A(^) = A+_ - A_+ (22) 

and the Heisenberg's exchange operator 

Ph = \{1 + (t^-ct2){1 + P,- P2)Pm (23) 

anticommute, where the operator exchanges the momenta(or coordinates). A*^^) 
violates the symmetry since the remaining part of the Hamiltonian and Ph commute. 

^ We have not assumed the counter correction which compensates the modification of the two- 
body propagator, since it represents some real effect of the crossed diagram. In QED for which the 
perturbative corrections work very well, we have to add the counter corrections to the vertex part. 
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For equal-mass limit, we have two equations. One is the equation (j^) with M = m 
and another is obtained by assigning a minus sign in front of A^^\ which is the equal- 
mass limit of the equation with m > M. It is the conjugate equation of (^) in the 
sense that @ is converted into it by the transformation Ph. It is easy to show that 
these equations have the common eigenvalue spectrum: If an eigenfunction Xn of ® 
belongs to some eigenvalue En, PhXu is the solution of the conjugate equation with 
the same eigenvalue. However, Xn does not have the definite P^-parity. 

For equal masses, it is reasonable to use the Salpeter equation which includes the 
projection factor (|^). The eigenvalue of this equation is different from the above En- 
The difference is, however, small since it is of the 4th order in the symmetry-breaking 
part of the Hamiltonian. 

4 Summary 

By assuming a proper boundary condition for the two-body propagation in the 
negative-energy state, we have proposed a new bound-state equation for the unequal- 
mass constituents. It has the symmetrical energy eigenvalues E and —E and reduces 
to the (one-body) Dirac equation in the infinite limit of one of the constituents masses. 
Secondly, we have discussed the normalization of the wave function and pointed out 
that the positive-definite probability density should be assumed. We can consistently 
calculate the observables of the bound state by assuming this normalization. 

The interaction Hamiltonian of the equation is constructed by diagonalizing the 
field theoretic Hamiltonian in the generalized Fock subspace of the two particles. 
Relativistic effects such as the retardation are taken into account in systematic way 
in the framework of the perturbation theory. 

We have further investigated the WFO's of the proposed equation in some detail 
for the instantaneous interaction and shown that the average WFO in the state 
which determines the leptonic decay rate of a pseudo-scalar meson is finite, which is 
in accord with the expectation from the field theory. 

One of the system for which we may utilize the two-body Dirac equation is the 
heavy flavored quarkonium. There are various possibilities of choosing the framework 
of the approximation according to the treatment for the high-momentum interaction. 
We have briefly discussed the correction factor to the leptonic decay width caused 
by the high momenta. The equation affords a good foundation of the heavy quark 
effective theory. 
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Appendix 

We examine the asymptotic {p oo) behavior of the momentum-space wave function 
and show that the average WFO Tr{757oV'(0)}/y2 is finite.0 

There are 4 partial amphtudes h^riip) in the ^5*0 state, with which the wave 
function is expanded as 

X(P) =EI:C.<(-P)V(P)'>.„(P)(^)'^ (24) 

where ci/2 = — c_i/2 = 1/v^ and the spinors u and f, for the particle 1 and 2 
respectively, are defined in 0. 

The average WFO for the annihilation decay through the axial- vector current is 
given by 

_LTV{,,,„,(0,}^-i-/(J^)./^ 



Y.{J{E^ + eM){E2 + r]m) - er]J{E^ - eM){E2 - vm)}K^ip)p^dp. (25) 



We first assume the Coulomb potential. The partial- wave equation for the ^5*0 
state is given by 



e'ri' 

p^E,{p)E,{q)E2{p)E2{qy 

J aI /1 2 I I aI 4 2 



{E - eE,ip) - 7]E2ip)}h,r,ip) = -e^ J dq 

e'-q' 

+ {ee'A\,^,,A'^^, + iqiAl,,A%_^^,}Qi{z)]Kj^,{q), (26) 

where z = {p^ + (f )/'^P(l and Qe{z) is the Legendre's function. Al^, is defined by 

Al,, = ,J{Ei{p) + em,){E,{q) + e'TTii). 

The asymptotic behavior of the wave function is determined from the integral 
region near the infinity. We then expand the both sides of ( ^6l) into the series of 
1/p and 1/q. We assume the power behavior of the wave function for large p. The 
independent amplitudes are chosen to be /^^(p) = J2ehee{p), hsip) = J^e^^eeip): 
hcip) = J2e^e-e{.p), and hoip) = Y^e^^s-sip)^ which are expanded, in the high- 
momentum region, in power series of 1/p: 



hx{p)=Y.Clp 



n ~l3x-1n-l 



n 



^See Rcf.|T^ and references therein, for the Salpeter equation. 
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Integrals on the right-hand side can be done if we neglect infrared- divergent terms 
which are irrelevant to the leading asymptotic behavior. Now, we can determine the 
asymptotic indices /?x's from consistency H: We get, for Ha and /i^, 



and 



2C\p-^- - EClp-P--' = -Cl-^ cot(^/?^)p-^- (27) 

TT I- Pa 2 

2C>-^- - EClp-^^-' = -C'.pjf^^ tM^B)?-'^, (28) 

TT [3b{2-Pb) 2 

where the terms of the higher power in 1/p are neglected. If we neglect the second 
term in the left-hand sides of (pT]), we find Pa in the range 1 < /^^^ < 2 and get 
Pb = Pa + ^ from (28). We obtain another series by neglecting the second term in 



For this, Pb is found to be in the range 2 < /3o < /3b < 3, where the lower 
bound Po corresponds to the upper bound i/n of a above which the index Pa from 
( p?!) becomes complex. Pa of the second series is given hj Pa = Pb + ^■ 

The asymptotic amplitudes he and ho are determined dependently on Ha and 
hs- We get, for the minimum indices 

Pc = miniPA + 2,Pb + 1) (29) 
Pd=Pa + 1. (30) 
We see that the average WFO ( p5l) is finite, because 

/3b > 1 and Pc > 2 

hold for the asymptotic amplitudes. This conclusion is valid even if the instantaneous 
exchange (transverse part) of the gauge boson is added. 
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